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Abstract 

We demonstrate a factorization formula for semi-inclusive deep-inelastic scattering with hadrons 
in the current fragmentation region detected at low transverse momentum. To facilitate the factor- 
ization, we introduce the transverse-momentum dependent parton distributions and fragmentation 
functions with gauge links slightly off the light-cone, and with soft-gluon radiations subtracted. 
We verify the factorization to one-loop order in perturbative quantum chromodynamics and argue 
that it is valid to all orders in perturbation theory. 
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I. INTRODUCTION 



In recent years, semi-inclusive deep-inelastic (SIDIS) lepton-nucleon scattering has 
emerged as an important tool to learn various aspects of perturbative and non-perturbative 
quantum chromodynamics (pQCD), the internal structure of the nucleon, in particular. The 
European Muon Collaboration experiment at CERN has provided us valuable information 
about the flavor dependence of quark fragmentation functions 0]. The HI and ZEUS col- 
laborations at the DESY HERA collider have measured the topology of the hadron final 
states in great detail and have compared them with the predictions of perturbative QCD0. 
In the area of polarized semi- inclusive DIS, the Spin Muon Collaboration, and recently the 
HERMES collaboration at DESY, have extracted the sea quark distributions and the polar- 
ized gluon distribution with controlled accuracy |3j, \$ . More recently, the target single-spin 
asymmetry measured by HERMES in semi-inclusive DIS is a new observable sensitive, for 
example, to the quark transversity distribution through the transverse momentum depen- 
dence of the produced hadron 0. 

In the semi- inclusive production of DIS, both the longitudinal momentum fraction z 
and the transverse-momentum Ph± of the hadron yield can be measured. When the trans- 
verse momentum is integrated over or when it is comparable to the hard photon-mass scale, 
Ph± ~ Qi the cross sections can be calculated from the standard pQCD formalism simi- 
lar to inclusive DIS and Feynman parton distributions jfj. In these cases, the theoretical 
tool has been well tested against experimental data with notable successes. When the 
transverse- momentum is much smaller than Q, but is still hard, Ph± 3> Aqcd, the cross 
section can be calculated again with integrated parton distributions augmented by small 
non-perturbative QCD corrections. The hard part contains the large double logarithms 
of the type a s In 2 Ph±/Q 2 - To make reliable predictions, these large logarithms must be 
summed 0, H| . An adequate formalism was developed by Collins and Soper in the case of 
e + e~ annihilation and shortly thereafter applied to the Drell-Yan process by Collins, 
Soper and Sterman (CSS) |l0|| . A first application of the CSS approach to SIDIS was made 
by Meng, Olness, and Soper Recently, a quantitative comparison between this theory 
and data from HERA collider has been made by Nadolsky, Stump and Yuan ^ . 



In this paper, we are interested in a special kinematic regime in SIDIS where Ph± is 
soft, i.e. on the order of Aqcd, and Q 2 is not too large, for example, on the order of tens 
or hundreds of GeV 2 . When Q 2 is large, the soft gluon radiations become important and 
can easily generate a large transverse momentum ^> Aqcd- Then the cross section for the 
hadron yield with Ph± ~ Aqcd is exponentially suppressed. To have a significant fraction 
of events with Ph± ~ Aqcd, fixed-target experiments with lepton beam energies on the 
order of tens to hundreds of GeV are preferred. The above kinematic regime is in fact ideal 
for studying transverse-momentum dependent (TMD) parton distributions in the nucleon 
and the related quark fragmentation functions. Recent interest in this subject has been 
stimulated by Collins's observation that semi-inclusive DIS at low-Pj. provides a tool to 
measure the quark transversity distribution jTj|]. The physics potential has been reinforced 



by the rediscovery of Siver's effect [1^ by Brodsky, Hwang, and Schimdt [1 



The main result of this paper is a QCD factorization theorem for the SIDIS cross section 
in the above kinematics region, accurate up to the power corrections {Ph±/Q 2 ) n and to all 
orders in perturbation theory. This factorization has been conjectured by Collins (Eq. 
(13)), following the early work of Collins and Soper on e + e~ annihilation 0. However, an 
exact statement of the factorization theorem requires an adequate definition of the TMD 
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parton distributions and fragmentation functions in QCD and a systematic factorization 
(and subtraction) of soft, collinear, and hard gluon contributions. In light of the recent 
development in this area (lR 17, 0, here we provide a first detailed examination of 



QCD radiative corrections in SIDIS, following the methodology of Ref. 

The factorization theorem we propose for the leading spin-independent structure function 

is 

F{x B ,z h ,P h± ,Q 2 ) = e l J d 2 k±d 2 p ± d 2 £ ± 

q=u,d,s,... 

xq (x B , k ± ,p 2 ,x B (,p) q T (z h , p ± , p 2 , (/ z h , p) S(i ± ,p 2 ,p) 
xH(Q 2 ,p 2 ,p)5 2 (z h k ± +p ± + £ ± -P h± ) , (1) 

where p is a renormalization (and collinear factorization) scale; p is a gluon rapidity cut-off 
parameter; the p and p dependence cancels among various factors. In a special system of 
coordinates in which x B ( = C/^h, one has C 2x B = s I z \ = Q 2 P- The physical interpretation 
of the factors are as follows: q is TMD quark distribution function depending on, among 
others, the Bjorken x B ; q is the TMD quark fragmentation function depending on, among 
others, the hadron momentum fraction z^, H represents the contribution of parton hard 
scattering and is a perturbation series in a s ; and, finally, the soft factor S comes from soft 
gluon radiations and is defined by a matrix element of Wilson lines in QCD vacuum. The 
above result shows that the hadron transverse momentum is generated from the combined 
effects of transverse-momentum of the quarks in the nucleon, soft gluon radiation, and the 
transverse-momentum of the quark fragmentation. There is no direct contribution from the 
TMD gluon distribution in this kinematic region. 

The main steps to establish the above factorization is as follows. In Sec. II, we introduce 
the TMD parton distribution and fragmentation function, and calculate them to one-loop 
order in perturbative QCD. The result contains both soft and collinear divergences and obeys 
the Collins and Soper evolution equation in the rapidity cut-off. We study the factorization 
of the TMD distributions by subtracting away the soft contributions. In Sec. Ill, one-loop 
result for semi-inclusive DIS scattering is obtained, and the factorization is shown to be true 
on the diagram-by-diagram basis. In Sec. IV, we generalize the one-loop result to all orders 
by identifying the leading regions for an arbitrary Feynman diagram using soft and collinear 
power counting. We then argue that a systematic factorization of the leading region leads 
to the general formula in Eq. (1). In Sec. V, the large logarithms in the perturbative 
expression are summed through solving evolution equations. We conclude the paper in Sec. 
VI. 

The factorization considered here can also be studied in the framework of soft-collinear 
effective theory developed recently in Refs. 0, 21, 2^, HI- We will leave this subject for a 
future publication. 



II. TRANSVERSE-MOMENTUM DEPENDENT PARTON DISTRIBUTION 

In the factorization formula (Eq. (1)), there is a factor q(x, k±, p 2 ,x B (, p) representing a 
transverse-momentum dependent (TMD) parton distribution, which differs from the usual 
Feynman parton distribution where the parton transverse-momentum has already been in- 
tegrated over. This object was introduced by Collins and Soper in the axial gauge and has a 
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number of interesting properties [9j|. In particular, it has a light-cone singularity and hence 
depends on the energy of the parent nucleon (related to Q i n addition to the parton's longi- 
tudinal momentum fraction x and transverse momentum k±. The sensitivity to the small-x 
gluon physics is controlled by a parameter p. 

In this paper, we follow a definition of TMD distribution in Feynman gauge with explicit 
gauge links jl6l |. We avoid the axial gauge because of the potential existence of gauge links 
at space-time infinity [l8j|. We calculate the TMDPD at one- loop order and show that 
it obeys the Collins-Soper evolution equation. The simplest definition of the distribution 
contains the soft-gluon effect which must be subtracted; we show how this can be done at 
the one-loop level. We also discuss theoretical difficulty to recover the integrated parton 
distribution from a direct transverse-momentum integration. 



A. Definition of TMDPD 



Consider a hadron, a nucleon for example, with four-momentum P. For convenience, 
we choose P along the ^-direction, P^ = (P°,0,0,P 3 ). In the limit P 3 — ► oo, the P^ 
is proportional to the light-cone vector (1,0,0,1). From now on, we use the light-cone 
coordinates k^ = (k° ± /c 3 )/a/2, and write any four- vector fc M in the form of (k~,k) = 
(k~ , k + , k±), where kj_ represents two perpendicular components (k x ,k y ). Let (xP + ,k_ s _) 
represent the momentum of a parton (quark or gluon) in the hadron. Let us start with the 
following definition of the transverse-momentum dependent quark distribution in a class of 
non-singular gauges 3, lit ]. 



Q(x, k±,n,x() 



2tt 



d 2 b 



(2tt 



-L e ib ± -k 1 



X (P 



V> g (r , 0, 6l)4(oo; r , 0, b ± h + £ v (w, 0)^(0) 



p 



(2) 



where the quark color indices are implicit, ip q is the quark field, v M is a time-like dimensionless 
(v 2 > 0) four- vector with zero transverse components (v ~, v + , 0), and C v is a gauge link along 



C v (oo; £) = exp ^— ig J d\v ■ A(Xv + . 



(3) 



The sign convention for the gauge coupling is = + igN 1 . As mentioned before, /i is 
an ultraviolet (UV) renormalization (or cut-off) scale. 

It is convenient to introduce the "gauge-invariant" quark field, 



The quark distribution becomes simply 



Q(x,k±,fi,xQ = - 



2tt 



-ix£-P + 



d 2 b A 
(27T) 



P 



^(r,o,&±)7 + ^(o) 



p 



(4) 



(5) 



The variable ( 2 denotes the combination (2P ■ v) 2 /v 2 = ( 2 . 

Physically, a parton interpretation of Q(x, k±, fi, () is the most natural if v is chosen along 
the conjugating light-cone direction of P M , i.e., ~ (1, 0, 0j_). However, as has been known 
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in the literature and re-emphasized recently by Collins jl9( , the distribution in this limit has 
logarithmic divergences (also called light-cone singularity) corresponding to contributions of 
virtual gluons with zero plus momentum £ + , or infinitely negative rapidity, \n£ + /£~. [In a 
physical process, the plus momentum of a parton is limited by the kinematics of scattering.] 
To avoid the divergence, we must introduce a rapidity cut-off for the gluons. One way to 
accomplish this is to introduce a non-light-like v^, such as with a d + ^ [16]. Then the 
contribution of the virtual gluons with rapidity smaller than \nv + /v~ is excluded from the 
parton distribution. As a consequence, a dimensional scalar ( 2 = (2P -v) 2 jv 2 emerges in the 
distribution. The limit of lifting the cut-off, v + — > 0, corresponds to £ — > oo. [In the following 
expressions, we will take this limit whenever we can.] The light-cone divergences are now 
reflected in the large logarithms involving (. The ^-evolution of the TMD distribution can 
be viewed as either the evolution in the gluon rapidity cut-off through v M or that in the 
energy of the incominghadron. The evolution is calculable in perturbation theory when kj_ 
is hard, i.e., > Aqcd M- 

Unless stated otherwise, we work in non-singular gauges, such as covariant gauges (in- 
cluding the Feynman gauge), for which the gauge potential vanishes at space-time infinity. 
In light-cone gauge, however, it is well known that the gauge potential is finite at infinity. 
Otherwise, the single spin asymmetry discussed by Brodsky et al. would disappear in such 
a gauge 15]. Belitsky, Ji and Yuan have shown that in singular gauges, one generally has to 
include gauge links at infinity 18]. These extra links can be found by imposing the gauge 
invariance of the parton densities starting from their definition in non-singular gauges. 

Since the two quark fields in Eq. (J2J) are separated along the spatial directions, the only 
ultraviolet divergence in Q(x, k±, /i, C) comes from the wave function renormalization of the 
quark fields and the gauge links. In this paper, we use dimensional regularization (DR) 
and modified minimal subtraction (MS) to treat ultraviolet divergences. If we use Eq. (J2J) 
naively in the axial gauge v ■ A — 0, then the ultraviolet divergence of Q(x, k±) is the same 
as the quark wave function renormalization in that gauge. Then the renormalization group 
equation becomes simple, 

V j- = 21fQ{x, k±, fi, xQ , (6) 

where is the anomalous dimension of the quark field in the axial gauge: 7^ = 
(3a s /4n)C F + (D(a 2 ). 



B. One-Loop Calculation 

In this subsection, we present the one-loop result for a quark TMD distribution in an "on- 
shell" quark. The calculation is important for a number of reasons. First, it shows clearly 
that the TMD distribution contains double logarithms in ( because of the collinear and light- 
cone divergences. It also serves as an explicit check for the evolution equation in rapidity 
cut-off. More importantly, the one- loop result allows one to devise QCD factorizations both 
for the distribution itself and for one-loop DIS cross section to be presented in the next 
section. 

We use a non-zero gluon mass A as an infrared regulator since there is no non-linear gluon 
coupling at one loop. One can use dimensional regularization beyond the leading order. The 
factorization is, of course, independent of the infrared regulator. Collinear singularities are 
regulated by none-zero quark masses. 
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(a) (b) (c) 

FIG. 1: Virtual- gluon contribution to one-loop transverse-momentum dependent quark distribution 
in an on-shell quark. The asymmetric diagrams from left-right reflection are not shown, but are 
included in the result. 



Let us first consider the virtual contribution shown in Fig. 1. For the self-energy diagram 
on the incoming quark leg, one has a contribution Q(x, k±) = 5(x — l)5 2 (k±)(Zp — 1) with 



1 + 



s^F 
An 



In4 + 2k4 



A 



(7) 



where m and A are the masses of the quark and gluon, respectively, Cp = (N 2 — 1)/{2N C ) 
with N c = 3; a term linear in N e = 2/e— 7^+ln 47r, where e = 4— d and 'jp the Euler constant, 
has been removed according to the MS scheme. The on-shell renormalization introduces the 
soft divergence in Zp, reflecting in the gluon mass dependence. For the self-energy on the 
gauge link, one has a similar contribution with Zp replaced by, 



CfLXs 
47T 



2 >4 



(8) 



Finally, the diagram with the virtual gluon vertex again has a similar contribution with Zp 
replaced by 



1 + 



a,C 



s^F 

An 



u 2 C 2 C 2 m 2 C 2 2tt 2 

21nA T + 21n^-ln 2 ^ T -21n^ r ln^- — + 4 



A 2 



(9) 



where we have made the approximation that 



A(P-v) 2 



is much larger than any other 



soft QCD scales. When £ is large, the double logarithms slow down the convergence of the 
pQCD series and call for a resummation which can be accomplished with the Collins-Soper 
equation (see the next subsection). In summary, the virtual diagrams give, 



Q(a?,fc_L,//,xC)|fig.i = S ( x ~ l)S 2 (k±){Z F + Z w + Z v - 3) , 



(10) 



where the dependence on soft scales m and A is implicit on the left. 

Now turn to the real gluon emission contributions shown in Fig. 2. The contribution 
from the diagram 2a) without the light-cone link, 



Q(X, fc±,/i,xC)|fig.2a 



a s Cp\l 



x 



2tt 2 



1 



k\ + x\ 2 + (1 
2xm 2 



x) 2 m 2 



(k\ + x\ 2 + (1 -x) 2 m 2 ) 2 _ 



(11) 
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(a) 



(b) 



(c) 



FIG. 2: Same as Fig. 1: real-gluon contribution. 

where we have taken e — ► 0. This must be done if we treat the TMDPD in the factorization 
formula as a physical observable. This, however, introduces certain problems in integrating 
out k± in DR, and we will discuss this more thoroughly in Sec. II. F. The transverse- 
momentum k± can go to zero, and therefore we cannot set the quark and gluon masses to 
zero too soon. However, the non-perturbative QCD physics will erase this sensitivity after 
factorization is formulated. 

The contribution from the diagram 2b), including the hermitian conjugation term, is, 



Q(x,fc_L,/i,xC)|fig.2b 



a s C F 



x 



7T" 



X 



1 



k\ + x\ 2 + (1 - x) 2 m 2 k\ + A 2 + C 2 (l - x) 2 



(12) 

The second term regularizes the light-cone singularity at x — 1. If one uses the usual 
regularization method of a plus function 24] and take the limit that ( is large, the above 
can be transformed into the following form, 



Q(X, fc±,/i,xC)|fig.2b 



a s Cf 



x 



71- 



(1 - x)+ k\ + x\ 2 + (1 - x) 2 m 2 



i oi s C F 1 
H — ^^-d(x — 1)— —m 



c 2 



2tt 2 



k 2 : + A 2 k 2 + A 



(13) 



which has a delta function at x = 1. 

Finally the contribution from the diagram 2c) with two gauge links is 



Q(X, fc_L,/i,xC)|fig.2c 



2|1 - x|C 2 



2vr 2 (A; 2 + A 2 + (1 - x) 2 C 2 )' 



(14) 



Formally, it vanishes when ( is large, except in the region where (1— x) 2 ( 2 is small. Therefore 
in the limit ( — ► oo, the above is the same as a 5 function at x — 1, 



Q(x, fc±,//,xC)|flg 



2c 



-<J(x 



2tt 2 



kj +A 2 



(15) 



We caution the reader that taking ( — >• cxd limit conflicts with the — > oo limit. In fact, 
in the above example, it turns a /c^-convergent integral into a divergent one. 



C. Collins-Soper Evolution in Hadron Energy or Gluon Rapidity Cut-Off 

As we have seen from the previous subsection, unlike the Feynman parton distributions 
which contain just the collinear singularities from the quark masses, the TMD distributions 
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contain in addition the light-cone singularities which are regulated by ( 2 . Since ( — > oo 
corresponds to the high-energy limit, ( dependence of the parton distribution is calculable 
in perturbative QCD, just like the renormalization scale dependence in /i. It turns out this 
is true only for large k±. 

The ^-evolution equation for Q(x, k±, /z, x() has been derived by Collins and Soper in 
the large ( limit jjj. Normally, we use Q(x, k±) in the small k± region. Let us extend this 
dependence to large k± and introduce the Fourier (or impact-parameter) representation, 

Q(x,b ± ,fi,x() = J <Pk ± e^ n± Q(x,k ± ,fi,x() . (16) 

The above integral should be convergent for non-zero b±. [There are UV divergences when 
6 = which we will not consider here.] The Collins-Soper evolution equation is 

C Q(x, b, n, xQ = {K{ii, b) + G(fi, xQ) Q(x, b, /I, x() , (17) 

where K depends on the UV renormalization scale /j, and infrared impact parameter b, and 
is non-perturbative when b is large; G is perturbative because /i and ( are hard; and both 
are free of gluon and quark mass singularity. The sum K + G is independent of UV scale /x 
and hence, 

li-fK = - lK = -li-^G, (18) 



where 7^ is the cusp anomalous dimension [25] and is a series in a s free of infrared sin- 
gularities. The derivation of the above equation in Feynman gauge has been given in Ref. 
jla lifj ]. In the above equation, any power correction of (Aqcd/C)" has been ignored and 
hence it is true only when ( ^> Aq CD . 

According to the result in the previous section, G gets a contribution from diagram lc 
only, whereas K gets a contribution from 2b. The sum is 

K(b, fi) + G(x(, 11) = In — 5— , (19) 

7T 4 

which is valid when b 2 is small and where 7^ is the Euler constant. The one-loop anomalous 
dimension is then, 

rv 

lK = ^2C F , (20) 

7T 

which is well known. Using the above renormalization group equation IjlSj) . one can sum 
over large logarithms ln( 2 b 2 in K + G when b is small (otherwise K is non-perturbative). 
Substituting the result into Eq. (fT7|). one finds a re-summed double-leading logarithms in 
(b (see Sec.V). 



D. Factorization of Soft Gluons in the TMD Parton Distribution 

From the viewpoint of QCD factorization, parton distributions are introduced to absorb 
collinear divergences when the quark masses are zero. From the one-loop result, it is seen that 
the TMD distribution contains both collinear and soft contributions, the latter is manifest 
through the gluon mass dependence. In this subsection, we will attempt to isolate and 
subtract the soft contribution from the above definition of TMDPD. 
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Let us first consider the self-energy diagram in Fig. 1. The diagram la contains the 
soft-gluon contribution which is obtained by making the soft approximation to the quark 
propagator and quark gluon coupling. Briefly, the soft approximation corresponds to ne- 
glecting the soft-gluon momentum in the numerator and the gluon momentum squared in 
the denominator; more discussion on the soft approximation is provided in Sec. IV. In the 
soft region where all of the components of £^ are small, the self-energy becomes 

dH l l 

(21) 



(4tt) 2 (£-p + ie) 2 (£ 2 - X 2 + ie) ' 

We can factorize the above contribution out of the parton distribution by subtracting it from 
the one-loop result of Zp. To make the subtraction mathematically convenient, we extend 
the soft approximation to the whole integration region of £ and use DR and MS scheme to 
get rid of the UV contribution. [There is no contribution from the collinear region because 
the integral is convergent in the massless quark limit.] The result is a contribution similar 
to the self-energy of an eikonal line. 

We can do the similar subtraction for diagram lb, by forming a soft approximation for 
the gluon interacting with the quark line 

_ 2 I ^_ 

19 J (2n) 4 (v£ + ie)(p-£ + ie)(£ 2 -\ 2 + ie) ' 1 ' 

with £ M restricted to the soft region. This time, however, the situation is more complicated. 
If one extends the integration £>* to all regions, there is also a collinear contribution coming 
from virtual gluons with momentum parallel to p M , as signified by the divergence of the zero 
quark mass. In other words, the simplified approach of subtracting away the whole integral 
will also take away a part of the collinear contribution. 

One may get around this by excluding the collinear gluon contribution with small £~ . 
This can be achieved by introducing a four-vector which has a large v + but relatively 
small v~ , and approximating the above soft contribution with the following integral, 

_. 2 f dH 2v -v 

~ %9 J (2tt) 4 (v ■ £ + ie)(v ■ £ + %e){£ 2 - A 2 + ie) ' ( } 

where we have replaced p by v. The above soft contribution includes soft gluons with £ + /£" 
limited by v + /v~ and t~ /£ + by v~ /v + . 

Subtracting the above contribution from the diagram lb, the remainder has a soft di- 
vergence in the gluon mass. This indicates that the soft and collinear divergences cannot 
be completely separated, as there are regions of loop momentum where collinear and soft 
divergences overlap. Therefore, one could in principle define Eq. (}2*2*j) with un-restricted 
^-integration as the "soft contribution". With this approach, the subtracted lb has no soft 
divergence. However, as we have mentioned before, the soft contribution is then not entirely 
soft. Two different approaches may be considered as two different subtraction schemes. 
Here, we use the first one. 

If we follow the above procedure, one finds the complete soft contribution in terms of the 
matrix element of Wilson lines, 

S(b ± ,fi 2 ,p) = ^(0|4,(6 ± ,-oo)4,^oo;fo ± )A ife (oo;0)A- fci (0;-oo)|0) , (24) 
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+ 











(a) 



(b) 



(c) 




FIG. 3: Soft-gluon contribution to the TMD parton distribution at one-loop order. The double lines 
represent eikonal line. 



where p = yv v + /v + v . we have made the color indices explicit (i,j,k,l = 1,2,3). The 
subtracted parton distribution can be defined as 



q(x, k ± ,n,x(,p) 



AT 

2tt 



d 2 b± ^,. S| (P|^(r,0,M 7 + ^(0)|P> 



(2tt)= 



S(b ± ,p 2 ,p) 



(25) 



This definition differs from that of Collins [Fj| in that the soft contribution here has no 
dependence in £~. Moreover, from our one-loop calculation, it is not clear that the above 
distribution has a well-defined limit when £ — > oo or v M becomes light-like, v 2 = 0, as claimed 
in 0. 

Let us calculate the one-loop soft subtraction, shown in Fig. 3. First, the diagrams with 
self-energy on all four of the Wilson lines, 



A soft g(x, /c_L)|dia g .3a = -S(x - l)5 2 (k ± )2(Z w - 1) 



(26) 



half of which cancels the self-energy of the gauge link in Eq. pOjl . 

The vertex correction of the gauge links, for which there is a factor of 2 to account for 
two virtual vertices, is 



A soft g(x, /c_L)|dia g .3b = -S(x - l)S 2 (k ± )2 



1 



1 



v-lv-lP-X 2 



S(x-l)6 2 (k ± )-^C F \n 



2n 1 v 2 v 2 



4(t> • v) 2 , p 2 
~ ln P ' 



(27) 



where the coefficient is just the cusp anomalous dimension |2f 

Now we consider the soft contribution from the real emission. Again there are two types 
of diagrams. The first type is the one with two gluons emitted from v or v. The contribution 
to the parton distribution is 



A soft g(a;, k ± )\ diag . 3c = 5(x - 1)—~G F 



71 



2 P k 2 +A 2 



(28) 



which, when integrated over k±, cancels out the self-energy contribution. The second type 
is the interference from the real gluon emission of the v and v lines. Plugging in the soft 
factor, 

a., _ , My- v) 2 1 



diae. 



,, = -5{x-l)—C F \n- y2 „ 2 2 



(29) 



When integrated over k±, it cancels out the vertex corrections. 



10 



E. Soft Contribution 



According to the definition of the soft contribution in Eq. (|24j). the one-loop result for 
the soft factor can be constructed from the above soft subtraction contributions, 



S(k ± ,n,p) = 5 2 (k ± ) 



S 2 (k ± ) + 



A soit q(x, k A 



5(x-l) 



2tt 2 



V 2 V 2 



When Fourier transformed to 6-space, it becomes, 



2tt 



-(2 - In//) In 



n 2 h 2 



(30) 



(31) 



This, however, cannot be used for 6 = 0, for which the integration over k± must be regularized 
so that S0= 0) = 1. 



F. Integrating over Transverse Momentum in TMDPD 

One would expect that after integrating over transverse momentum, a TMD parton dis- 
tribution reduces to the usual Feynman parton distribution. This, in fact, is not straight- 
forward. 

First of all, we have chosen DR and MS to regulate ultraviolet divergences. In the 
TMDPD, the e — > limit and MS subtraction have already been performed as it represents 
an observable in 4-dimension. On the other hand, a Feynman parton distribution is ob- 
tained first by integrating the transverse-momentum and then performing UV subtraction. 
Since the integral is a divergent one, the procedure of integration and subtraction is not 
interchangeable. To get around this, one may use a momentum cut-off to regularize the UV 
divergence but this is hard to implement beyond one- loop without destroying the gauge 
symmetry. [A consistent regularization might be a discrete space-time lattice, however, in 
practice this is hard to implement in Minkowski space.] 

Even when there is a consistent cut-off regulator, one may still have a problem with the 
light-cone singularities. While the gauge link in Feynman distributions is along the light- 
cone, we have chosen an off-light-cone gauge link to regulate these divergences. Therefore, 
one cannot expect that after integrating over a TMDPD with a non-light-like gauge the 
Feynman distribution recovers. 

What happens if one takes the light-cone limit of the non-light-like gauge link after 
integrating over the transverse momentum in a TMDPD? The standard integrated parton 
distribution still does not emerge if dimensional regularization is used. If one chooses the 
non-light-like gauge from the beginning, the gluon propagator has a term proportional to 
v 2 (if v is the direction of the gauge link). This term contributes to the TMDPD in loop 
calculations. When the light-cone limit, v 2 — > 0, is taken, these contributions would vanish 
if there were no axial-gauge singularities at k ■ v = 0. In practice, however, the singularities 
are present and the limit v 2 — > does not reproduce the result obtained with v 2 = set in 
the beginning in dimensional regularization. 

To summarize, it is nontrivial to recover a Feynman parton distribution by integrating 
over the transverse-momentum in a TMDPD. One could cut off the integral by hand, but 
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a UV regularization scheme must be used which implements the same cut-off in the loop 
integrals. This is difficult to construct beyond one loop. For the same reason, the light-cone 
limit of the vector v M is not analytical. 



G. TMD Fragmentation Function 

The transverse-momentum-dependent fragmentation function (TMDFF) has a similar 
definition as the TMD parton distribution. Many of the results discussed in the previous 
subsections can be immediately translated into those for TMDFF. Here, we sketch the main 
results briefly. 

Use qh(z, kj_, p, Q as a notation for a subtracted quark fragmentation function into the 
hadron h, 



q h (z,P h± ,p,(/z,p) = 1 /" f^#V^--^) (32) 



Pi 

2z J 2tt (2tt) 2 ~ 
x E I D ^- 00 ' 0)^(0)^^)7^ 

X a 

x (P,X|(VUr , b)4(CX -oo)|0)/S(6j_, n, p) , 

where v is mainly along the light-cone direction conjugating to Ph', k + = P^/z and k± = 
—Ph±/ Z ', an d a is a color index. The variable ( is defined as 

C 2 = A(P h ■ vf/v 2 . (33) 

For a quark fragmenting into a quark, the leading order result is normalized to q(z,P±) = 
5(z-l)5 2 (P ± ). 

It is not difficult to see that the quark fragmentation function in a quark can be obtained 
by a simple substitution of the corresponding quark distribution in a quark, 

1 / 1 P 

q h (z,P ± ,fi 2 ,(/z,p) = -q[-,—,p 2 ,(/z,p 

z \z z 

= q(z,P ± ,p 2 ,z(,p) , (34) 

where the second equality holds at one-loop order. Moreover, q^ satisfies the same Collins- 
Soper equation in ( evolution as that of q in ( evolution. 



III. ONE-LOOP FACTORIZATION 

In this section, we show that the factorization formula, Eq. (1), is valid at one-loop order. 
To accomplish this, semi-inclusive DIS on a single quark target is studied. The result can 
be easily translated into that for a non-perturbative hadronic target. 

In the first subsection, we establish notation and normalization for the tree level result. 
In the second subsection, we state and explain the content of the factorization theorem. In 
the following two subsections, we will verify its correctness for a single quark target on the 
diagram-by-diagram basis: first for the virtual corrections, and then for the real corrections. 
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A. Notation and Tree Normalization 



We choose a coordinate system for semi-inclusive DIS in which the nucleon is travel- 
ling along the ^-direction. Introduce the light-cone vectors (p° , p x , p v , p z ) = A(l, 0,0,1), 
(n° , n x , n v , n z ) = (1, 0, 0, — 1)/2A, and p ■ n — 1, where A is an arbitrary parameter. The 
initial nucleon momentum P can be written as, 

P» = p»+ (M 2 /2)n fl , (35) 

where M is the nucleon mass. The photon momentum is q = £ — £', where £ and £' are the 
initial and final lepton momenta, respectively. We choose the photon momentum in the 
negative- z direction, 

= + ^U" , (36) 

where £ ~ xb = Q 2 /2(P ■ q) when M 2 /Q 2 is neglected and Q 2 = —q 2 . The so-called hadron 
frame is obtained by making a particular choice of A 0. 

One has the option of either fixing the lepton plane as the xz plane or the hadron plane 
as the xz plane. In either case, there is an azimuthal angle between the two planes, and 
for the simplicity of our discussion, we integrate out this angle. The detected hadron has 
a momentum P^ mainly along the n 11 (negative z) direction with z^ fraction of the photon 
momentum component in the same direction, and with transverse momentum Ph± which is 
invariant under the boost along the z-direction. As indicated earlier, Ph± is considered to 
be soft (on the order of Aqcd)- If Ph± Aqcd, a different factorization formula exists in 
which only the integrated parton distributions and fragmentation functions enter. 

The semi-inclusive DIS cross section under the one-photon exchange is 

- 2naL y£^(P, q ,P h ), (37) 



dxBdydzhd 2 Ph± Q 
where the unpolarized lepton tensor is 

t v = 2(ff + ff -g^Q 2 /2) 

= (Q7^)(l- y + y 72)HVn + --- , (38) 

where y is the fraction of the lepton energy loss, 1 — E'/E. Since we are going to integrate 
over the azimuthal angle 0, only the structure = g^ u — p fl n l/ — p u n fl will survive. 
The hadron tensor has the following expression in QCD, 

W^(P,q,P h ) = T^E/ ^e ig< (P\UO\XP h }(XP h \J u (0)\P) , (39) 

where J M is the electromagnetic current of the quarks, X represents all other final-state 
hadrons other than the observed particle h. The variable Zh can be defined as P ■ Ph/P ■ q 
ovP'/q-. 

A simple calculation on the single-quark target yields that, 

= -\g^8{x B - \)8{z h - l)5 2 (P h± ) + ... . (40) 
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In this case, it is known 

q (0 \x B ,k ± ) = 5(x B -l)5 2 (k ± ) , 

q {0 \z h ,p ± ) = 8(z h -l)S 2 (p ± ) . (41) 
It is easy to translate the above into a result for a physical hadron, 

W> w = ~g^ J d 2 k ± q(x B ,k ± ) j d 2 p ± q(z h ,p ± )5 2 (z h k ± +p ± ~P h± ) . (42) 

Therefore the cross section is, 

x J d 2 k±q(x B ,k ± ) J d 2 p ± q h (z h ,p±)5 2 (z h k ± + p± - P hA _) , (43) 

where s = (P + £) 2 , and we have kept only the ^-independent term. This result is known in 
the literature 12811 . 



B. General Form of Factorization 

In the following discussion, we are interested in the leading structure F(x B , Zh, Ph±,Q 2 ) 
only, 

= -lfi;F{xB, z h , P h ±, Q 2 ) + ... . (44) 

The other structures factorize in a similar way. The form of the factorization theorem we 
want to show is 



F(x B ,z h} P h± ,Q 2 ) = ^2 e l / d 2 k±d 2 p±d 2 

j „ J 



q=u,d,s,. 



xq (x B ,/c_ L ,/i 2 ,x B C,p) Qh [z h ,p ±J fi 2 ,(/z h ,p) S{l x ,lJ?,p) 
xH(Q 2 ^ 2 ,p)5 2 (z h k ± +p ± +£ ± -P h± ) , (45) 

where in a special system of coordinates: ( 2 = (Q 2 /x 2 B )p and ( 2 = {Q 2 z 2 l )p and p = 
■\/v~v + /v + v~ is a gluon rapidity cut-off parameter. The above result is accurate up to 
powers in (P^ ± /Q 2 ) n for soft Ph± ~ Aqcd- There is no direct contribution from the gluon 
distribution in this kinematic region. There is no convolution involving the longitudinal 
momentum fractions, x B and Zh, typical in other hard processes. The transverse-momentum 
integrals show that the hadron transverse-momentum can be generated from the initial state 
parton, final state fragmentation, and the soft gluon radiation. 

The renormalization and collinear factorization scale p cancels among the four factors, as 
the structure function should be /^-independent. As in the inclusive case, one can choose p? = 
Q 2 to eliminate the large logarithms in the hard factor. The soft-factorization parameter p 
depends on the directions of the Wilson lines and must also be cancelled among all factors. 
The Collins-Soper equation allows studying the double logarithmic dependence of Q 2 in the 
distribution and fragmentation functions. 
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FIG. 4: One-loop virtual correction to semi-inclusive DIS. 



In principle, for Ph± ~ Q, the above factorization formula breaks down because of the 
power corrections, and because the transverse-momentum Phj_ is now mainly generated from 
multi-jets production. However, it is convenient to extrapolate the above factorization to 
all Ph±, and introduce the impact-parameter space representation, 



F(x B ,z h ,b,Q 2 ) = ^2 e g<? {xB,z h b,[i 2 ,x B (,p) q (z h ,b, p 2 ,(/z h , p^J 



q=u,d,s, 

xS(b,p 2 ,p)H(Q 2 ,p 2 ,p) 



(46) 



The convolution in the transverse momentum becomes a product of Fourier factors. 

At tree level, S^(b, p 2 , p) = 1, and = 1. Let us show that at one-loop, the above 
factorization is still valid, and calculate the correction to the hard part at this order. 



C. Virtual Corrections 

We first consider factorization of one-loop virtual corrections to the tree process. The 
factorization actually holds diagram by diagram, and therefore we will study the momentum 
flow in individual diagrams and extract the corresponding hard factor. 

Three diagrams shown in Fig. 4 correspond to the initial and final state wave function 
renormalization and vertex corrections. 

The self-energy correction is straightforward, 

F = 5(x B - \)8{z h - l)5 2 (P h± )(l + 2(Z F - 1)) + ... (47) 

where Zp is given in Eq. (JJJ). Zp contains both soft and collinear contribution in the on- 
shell scheme. If the self-energy is associated with the initial state quark, the collinear part 
of Zp is absorbed by the self-energy correction on the parton distribution, corresponding to 
diagram la subtracted by diagram 3a. The remaining soft contribution is attributed to the 
soft factor shown in diagram 3a. There is no contribution from the fragmentation function, 
nor is there a contribution from the hard part. A similar argument shows the self-energy 
correction to the final-state quark can be absorbed by the fragmentation function and the 
soft factor, yielding no contribution to the hard part. 

The vertex correction produces exactly the same expression for F as the above, except 
Zp is replaced by, 

Z r = 1 -^C F L91 + ^% + 2^^-4^-^] (48) 
Air \ p z m z A z mr m z 3 / 
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FIG. 5: One-loop real correction to semi-inclusive DIS. 

The UV divergence in the above expression cancels that in Zp because the sum of one-loop 
virtual corrections has no UV divergence. The above result contains a quark distribution 
part shown in Fig. lc subtracted by Fig. 3b, and a fragmentation function with a similar 
structure, and a soft contribution in Fig. 3b. Subtracting all of the above from Zy — 1, we 
find a left-over hard contribution, 



ol. 



H^(Q 2 ,p 2 ,p) = ^C F 

Z7T 



(1 + lnp 2 ) In 



Q i 2 i 1 i 2 2 i 2 

— - In p £ + - ln z /T + 7T 
p l 4 



(49) 



where we have chosen a coordinate system in which x B ( 
on the quasi-light-like vectors v and v is simply through a combination, p 
The dependence on Q 2 is of the form of single logarithms and can be controlled by a renor- 
malization group equation because it contains no additional scale other than p. 



(/zh and therefore the dependence 
sjv ~v + /v + v~. 



D. Real Corrections 

The one-loop real corrections are shown in Fig. 4. Since we are interested in the topol- 
ogy of the final state in which the struck quark carries the dominant part of the energy- 
momentum of the current region, the emitted gluons are considered to be either soft or 
collinear. Therefore, there is no contribution to the hard scattering kernel from any of these 
diagrams. Our job is to show that these diagrams can be properly taken into account by 
the known one-loop parton distribution, fragmentation function, and the soft factor. 

Let us start with the ladder diagram shown in Fig. 5a. The soft-gluon radiation generates 
a transverse-momentum for the struck quark. There is no contribution from the fragmenta- 
tion function because the contribution from the final state with a gluon in the n M direction 
and a soft quark is power suppressed. Therefore, the diagram must be factorizable into the 
parton distribution in Fig. 2a subtracted off Fig. 3c, and the soft factor in Fig. 3c. A simple 
calculation of the diagram yields, 



-±C F S(z h - l)(l-x B ) 
2xm 2 



1 



(P 2 ± +x B \ 2 + (l-x B ) 2 m 2 ) 2 



P 2 ± + x B X 2 + (l-x B ) 2 m 2 

(50) 



Indeed, the above expressions can easily be reproduced by the factorization formula with a 
one-loop result for and and tree level and H( \ 
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Similarly for Fig. 4b, we have 



2^ 



C F 5{x B -l){\ 



1 



P 2 ± + z h X 2 + (1 - Z h ) 2 m 2 



2xm 2 



(P 2 ± + z h X 2 + (l-z h ) 2 m 2 ) 2 



(51) 



which again can be reproduced by the factorization formula with the one-loop fragmentation 
function and the soft factor S, and the tree-level parton distribution and the hard part. 

Finally, let us consider the diagram Fig. 5c and its hermitian conjugate. In the region 
where Ph± is small, we find three distinct contributions: 



a s Ci 
~2~^ 



■6(z h -l)- 



2x B 



CU.Cf 



2tt 2 
a s C F 



${Xh 



1 - x B )- 
, 2z h 



1 



p2 



1 



Zh 



+ X B X 2 + (1 — x B ) 2 m 2 

1 



p2 
r h±_ 



2ir 2 



25{x B - l)8(z h - V 



PL + A 2 



z h X 2 



In 



1 - z h ) 2 m 2 



Q 2 



PL + a 2 



(52) 



where the first term corresponds to a gluon collinear to the initial quark, the second term a 
gluon collinear to the final state quark, and the third term a soft gluon. All these terms are 
reproduced by the factorization formula with one-loop parton distribution, fragmentation 
function, and the soft factor. 

Therefore we conclude that at the one-loop level, the general factorization formula holds. 



IV. FACTORIZATION TO ALL-ORDERS IN PERTURBATION THEORY 

In this section, we argue that the factorization formula we stated in the previous section 
holds to all orders in perturbative QCD. To make such arguments, we follow the steps 
outlined in an excellent review article by Collins, Sterman, and Soper j2^|. One must 
consider a general Feynman diagram and study its leading contributions to the SIDIS cross 
section. The contributions from different regions of the internal momentum integrations 
are characterized by the reduced diagrams which correspond to pinched surfaces in the 
space of integration variables. The leading reduced diagrams can be determined by infrared 
power counting. The remaining steps involve decoupling the Lorentz and color indices, and 
using soft approximation to disentangle momentum integrals in the different parts of the 
reduced diagrams. To simplify the derivation, one must use the (generalized) Ward identities 
extensively. After decoupling and replacing the various factors by the parton distribution, 
fragmentation function and the soft function, one finally arrives at the general form of 
factorization. 

Using the fact that the physical observables are independent of renormalization and soft- 
collinear factorization scales, large double and single logarithms involved in the factorization 
formula can be summed. The final expression is useful to describe experimental data when 
combined with result from perturbative calculations of a fixed order. 
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FIG. 6: A general reduced diagram for semi-inclusive DIS. 



A. Reduced Diagrams, Power Counting, and Leading Regions 

The contribution of an arbitrary (cut) Feynman diagram to the SIDIS cross section can be 
classified in term of pinched surfaces corresponding to the solutions of Landau equations [^U 
EuL IHf • Coleman and Norton observed that these pinched surfaces can be pictured in terms 
of physical space-time processes (or reduced diagrams) |33]. According to the kinematic 
constraints of SIDIS, it is not difficult to see that the most general reduced diagrams have the 
structure shown in Fig. 6, in which the initial nucleon evolves into a target fragmentation jet 
Jt plus a set of collinear quarks and gluons (solid lines) entering the hard-interaction vertex 
with the highly virtual photon 7*. A new set of collinear quarks and gluons (solid lines) 
emerges from the hard vertex in a new direction (opposite direction in the collinear frames), 
and fragments into the observed hadron and the unobserved jet J c . Fig. 6 is actually a cut 
diagram including the complex-conjugated amplitudes, corresponding to the measured cross 
section. Therefore, we will use additional indices L and R to label jets on the left and right 
sides of the cut (indicated by the vertical dashed line), respectively. For example, J c r labels 
the current jet on the right-hand-side of the cut. In addition, there is a soft sub-diagram S 
with soft quark and gluon lines (shown by dashed lines) connecting the jets and hard parts. 

Let us count the degree of infrared divergence 00(G) of each reduced diagrams G. It can 
be constructed from the sum of the degrees of divergences for the jets and the soft part, 



The power counting for the soft function is straightforward: If we use E b and E* to denote 
the number of soft boson (gluon) and fermion external lines, then it is well known that 



from a simple dimensional analysis in coordinate space. Note that us includes the propaga- 
tors of the external lines and the associated integration measure. 

Let us use pj to denote the number of collinear quark or gluon-with-physical-polarization 
lines entering the hard part from jet J; use lj to represent the number of collinear gluons 
of longitudinal polarization through a similar attachment; use E b j^ to denote the number 
of soft boson or fermion lines connecting the soft part to the jets; use EJ L r to label the 



u{G) = oj JtL + u JtR + uj cL + lo JcR + u s • 



(53) 




(54) 
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FIG. 7: The leading region for semi-inclusive DIS. 



number of soft bosons or fermions connecting to the left or right hard part; and finally use 

(3) 

v j to label the number of three-point vertices in the jet, and sj the number of soft gluons 
with scalar polarization attaching to the jet. Then it is easy to see that 



E b = E b JtL + E\ R + E\ L + E\ R + E h HL + E b HR , 

Ef = E L + E L + E L + E L + e hl + E L ■ (55) 

The soft power associated with each collinear jet is 

loj = 2Lj -Nj + tj , (56) 

where Lj is the number of loops in the jet (each contributing two powers), and Nj is the 
number of internal lines (each contributing one power), and tj is the numerator suppression 
factor which in Feynman gauge is equal to max[t> ^ — lj — sj,0]/2 j^J]. The number of loops 
can be calculated using 

Nj - vf - vf = Lj , (57) 

where v j 3,4 "* are the number of three-and four-point vertices, respectively. The relation 
between the number of vertices and lines is, 

3vf + Avf +pj + lj = 2Nj +E b J + E f J + Ij, (58) 

where lj is the number of (initial and final) external lines in the jet. From the above, it is 
easy to see that 

= ^pj-s J -E b J -E f J - lj) + l -[sj + lj- vf) 9 (sj + lj- vf) . (59) 
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Combining the results from four jets, one finds, 

u(G) > -(pj tL +Pj tR +Pj cL +Pj cR )-^(It + I c ) 
+ - (E b r + E b } + E b j + E b ) 

2 V J tL JtR JcL JcR) 

1 / 

~2 ( Sj " + Sj m + S Jcl + S JcR> 

+ E L + E{ R + Ei L + E{ r 
+E b HL + E b HR + \{E f HL + E f HR ) 

+\{sj tL + i JtL -vf L )e{s JtL + i JtL -vf L ) 
+\{sj tR + h tR -v^)e( SJtR + i JtR -vf tR ) 
+l{^ + h cL -v i Z)e{ SJcL + i JcL -vfl) 

+ \ (s JcR + l JcR - vf cR ) 9 (s JcR + l JcR - vfl) . (60) 

From the above, the largest possible degree of infrared divergence is if the initial and final 
state hadrons are replaced by a perturbative parton (I t — I c — 1). 

According to the above result for u>(G), leading reduced diagrams (leading region) must 
satisfy the following conditions: 

• No soft fermion lines, 

• No soft gluon lines attached to the hard parts, 

• Soft gluon lines attached to jets must be longitudinally polarized, 

• In each jet, one quark line plus an arbitrary number of longitudinally polarized gluons 
attached to the corresponding hard part, 

• The number of three-point vertices in a jet must be larger or equal to the number of 
soft and longitudinally polarized gluon attachments. 

In Fig. 7, we show the leading reduced diagrams satisfying the above conditions. As 
indicated already, the collinear gluons are longitudinally polarized. 



B. Factorization of Collinear Gluons 



Let us first factorize the longitudinally-polarized colliner gluons from the hard parts. 
This can be done using the approach discussed in 29j | . For definiteness, let us consider 



the collinear gluons from the initial state nucleon. Because the gluons are longitudinally 
polarized, the gluon gauge potential can be replaced by 

A» = A ■ np" , (61) 
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where p^ is the light-cone momentum to which the initial nucleon momentum is proportional. 
The effects of these gluons on the hard part can be factorized through the Ward identity, 

(f\Td^^ 1 )d^(^)...d IMl A^(Q\i) = , (62) 

where \i) and |/) are physical states. Applying this identity repeatedly leads to the conclu- 
sion that the collinear gluons can be viewed as attaching to an eikonal line in the conjugating 
light-cone direction n^. This result can be understood in an intuitive way: The longitudi- 
nally polarized gluons cannot resolve the internal dynamics of the hard scattering. It can, 
however, be sensitive to the overall flow of the color charge. The hard interaction is a light- 
cone dominated process along the n M direction in the coordinate space. This is also the 
direction along which the final state jet is formed. Thus the collinear gluons mainly scatter 
with the color-charge flow in this direction. 

The Feynman momentum a; of a collinear gluon has a lower limit in a physical process. 
For example, the smallest x that a gluon may have is on the order of M / Q. Only in the limit 
Q — ► oo, can there be near zero-momentum gluons participating in the scattering. When a 
collinear gluon has a small x, its light-cone energy is large, and its contribution to the cross 
section can be calculated perturbatively. Therefore, one can introduce a parameter that 
separates contributions of the gluons with different rapidities. The collinear gluons with x 
larger than a certain cut-off are included in the parton distributions; others are included in 
the hard factor. Of course, the physical cross section is independent of this parameter. In 
the inclusive case, the singular contribution from small-x gluons cancels between the real 
and virtual diagrams. 

To define a parton distribution with virtual gluons of limited rapidity, one can introduce 
a rapidity cut-off. The most straightforward approach is to implement a lower cut-off in x. 
A more convenient approach, as we discussed in the one-loop case, is to introduce a quasi- 
light-cone vector, v M , which is close but not exactly in the n M direction, and to assume that 
all collinear gluons couple to a colored jet moving in this direction. It can be checked that in 
this approach only the gluons with k + /k~ > v + jv~ are included in the parton distribution. 

The collinear gluons from J t can be factorized in a similar way. Here a quasi-light- 
cone vector v must be introduced to limit the small- a; gluon contribution to the TMD 
fragmentation function: only collinear gluons with k~ /k + > v~ /v + are included in the non- 
perturbative function. The left and right parts of the cut diagrams can be treated in a 
symmetric way. 

C. Soft Approximation and Soft Factor 

The soft gluons are attached to the target and current jets, and can be factorized using the 
Grammer-Yennie approximation j34| (or soft approximation). The approximation consists 
of two steps. The first step is to neglect any soft momentum in the numerators of the 
jet factors. One of the consequences is that the gluon polarization is effectively along the 
conjugating light-cone direction of the jet (longitudinally polarized). The second step is to 
neglect k 2 compared to k-nk-p in the jet denominator. This approximation is not uniformly 
true in the soft region. In fact, in the so-called Glauber region, where k\ k ■ nk ■ p, the 
approximation fails |35[ . If, however, the momentum k + or k~ is not trapped, one can deform 
the contour integration to a region where k ■ nk ■ p ^> k\, so that the approximation can 
still be applied. It is known that for semi-inclusive hadron production in e + e~ annihilation, 
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FIG. 8: The leading region for SIDIS after soft and collinear factorizations. 



the deformation can be easily performed |29j. In inclusive Drell-Yan, this happens only 
after summing the final state interaction diagrams j3(|. In the present case, the soft gluons 
interact with the current and target jets; all of these interactions are final state interactions. 
Hence, all physical poles appear in the upper-half plane. As such, the contour deformation 
can be done straightforwardly. 

After the soft approximation, one can again use the Ward identity to factorize all the soft 
gluons from the jets. The physical effect of a jet can be replaced by a Wilson line along the 
jet direction. Again to avoid the light-cone singularity, the Wilson line can be chosen to be 
off the light-cone along the v or v direction. After factorizing the gluons from the jets, and 
summing over all soft contributions, a soft factor emerges: 

S(b ± ,p 2 ,p) = ^Tt<0|4(6 ± ,-(X3)4(oo;6j.)£ w (oo;0)£ e (0;-oo)|0> , (63) 

which appears as a factor in the factorization theorem. Now the leading region has the form 
shown in Fig. 8. 

The soft factor is renormalization-scale dependent. The renormalization group equation 

is 

dS(b ± ,p 2 ,p) 



/'- 



dp 



7s(p)S(b±,p ,p) 



(64) 



where 75 (p) is the anomalous dimension of the Wilson lines in the definition. At one- loop 
order, one has 

7 5 = -^[2-lnp 2 ]+... (65) 



7T 



which is p dependent. The anomalous dimension at higher-order has been studied in Ref. 



D. Subtracted and Unsubtracted Parton Distributions and Fragmentation Func- 
tions 

Let us consider the target jet factor which has been factored from the hard part, with the 
soft factor factorized out as well. It is shown on the left-hand side in Fig. 9. The internal 
loop momenta of the gluons are restricted. Their k + components have a lower limit because 
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FIG. 9: Relation between soft-subtracted (left-hand side) and un-subtracted parton distributions. 
The subtracted distribution has a light-gray blob and the un-subtracted one has a dark blob. The 
denominator on the right-hand side is a soft factor. 



of the gauge- link direction v. The k~ components also have a lower limit: when the soft 
gluons are factored, the gluons with k~/k + smaller than v~/v + have been factored out of 
the jet. 

Therefore the jet factor is not the same as the parton distribution Q(x, b, p 2 , xQ defined 
in Sec. II. Rather, it is the same as the soft-subtracted parton distribution q(x,b, p 2 ,x(, p). 
The relationship of the two is shown in Fig. 9. 

The renormalization group equations for Q and Q are known. After subtracting the soft 
factor, the equation has to be modified by including the anomalous dimension for the soft 
factor: 

dq(x,b,p 2 ,x(,p) /\\ / i 2 /- \ (RR \ 

V ^ = (27f - 7s(p)J Q{x, b, p , xQ, p) . (66) 

A similar equation holds for the TMD fragmentation function. The Collins and Soper 
equation for q is the same as for Q: 

d 

C^O, b, p 2 , x(, p) = [K(n, b) + G(p, x())q(x, b, p 2 , x(, p) . (67) 
We will solve this equation to sum over large logarithms in Sec. V. 



E. Subtraction Method 

For a given Feynman diagram, there are multiple leading regions. To ensure the factor- 
ization works, one, in principle, has to supply a subtraction method which allows separating 
contributions from different leading regions. In particular, the subtraction method must 
provide a systematic way of handling the overlapping contribution of different leading re- 
gions. 

The easiest way to develop a subtraction method in gauge theory is to choose the axial 
gauge. For example, the factorization for inclusive DIS in the axial gauge can be developed 
using the Bethe-Salpeter formalism and has been used to calculate the anomalous dimension 
of parton distributions at two-loop order j3^j . For factorization involving collinear and soft 
divergences, a subtraction method in the axial gauge has been developed by Collins and 
Soper [^]. 

In covariant gauge, a systematic subtraction is complicated, and has not yet been fully 
developed in the literature. At the one-loop level, an example has been provided by Collins 
and Hautmann 17]. We have checked that the new subtraction method corresponds to a 
particular choice of p in this paper. It would be interesting to pursue this subtraction to 
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higher-order. In particular, higher-order calculations help to clarify the roles of different 
non-perturbative matrix elements which have the same one- loop result. 

The best approach to treating overlapping infrared divergences in a multi-loop case might 
be the soft-collinear effective theory mentioned in the introduction. Here we assume this 
can be done in principle and leave a more careful discussion for future publication. 



F. Factorization and p-Independence 

Collecting all factors in Fig. 8, one finally has the following factorization formula: 
F(x B ,z h ,b,Q 2 ) = ^2 e ?<? {xB,z h b,p 2 ,x B (,p) q (z h ,b, p 2 ,(/z h , p") 

q=u,d,s,... 

xS(b,p 2 ,p)H(Q 2 ,p 2 ,p) , 



(68) 



advertised earlier. 

It can be shown that the above expression is independent of p. At one-loop order, this is 
easy to see, 



dS(b, p) a s C F ( p 2 b 2 2 \ 
p - B — = -—ln\- r eT) S(b,p), 



On the other hand, 



dS(b,p)q(x,b,p) _ a s C F ^ f pQ 2 b 2 p27g _ 1 



dp 



2tt 



4 



S(b,p)q(x,b,p) 



If the structure function is independent of p, the above requires H(p) to evolve in p, 



dH(p) _ a s C F 
dp 7i 



In 



P& 
P 2 



- 1 



H(P) 



(69) 



(70) 



(71) 



It is easy to check that our one-loop hard part satisfies the above equation. 

At higher orders, the p independence is guaranteed because one can view the factorization 
formula as a definition for the hard part. 



V. SUMMING OVER LARGE LOGARITHMS 

From the factorization formula and the evolution equations, we can get an expression for 
the structure function in which the large logarithms involving momentum Q are summed 
over. 

First of all, there are large logarithms in K + G (which is independent of the renormal- 
ization scale). To sum it, we solve the renormalization group equation to get 

C m da 

K(b, p) + G(x(, p) = K(b, p L ) + G(x(, p H ) - / -4ikHp)) . (72) 

To isolate the large logarithms, one has to choose p L to be on the order of A QCD and p H to 
be on the order of (. Therefore, we let 

p L = dM N ; p H = C 2 x( = C 2 Q^p , (73) 



24 



where is the mass of the nucleoli. 

Substituting the above into the Collins-Soper equation for q(x,b, p 2 ,x(, p), the large 
logarithms in ( can be factorized, 



q(x,b,p,x(,p) = expl- 



p [ \ p J 

xq(x,b,p,x( = p L /C 2 ,p) , 



7x(«(/i)) - K(b, p L ) - G(p/ C 2 , p) 



(74) 



where the exponential factor contains the entire dependence on (, in particular, the large 
Sudakov double logarithms. However, the above expression contains much more than just 
the leading double logarithms; it contains all the sub-leading logs as well. 
Similarly, one can find the solution for the fragmentation function, 



q(z,b,p,(/z,p) = exp 



r c ^ z dp 



In 



C2C 
zp 



1k{&{p)) - K(b,p L ) - G(p/C 2 ,p) 



xq(z,b,p,( /z = p L /C 2 ,p) 



(75) 



If we choose a frame in which xb( = C/ Z h, then the exponential factors in q and q become 
the same, and moreover ( 2 x 2 B = C 2 1 z h — Q 2 P- 

Let us study the renormalization group equation for the hard part. The physical cross 
section is, of course, independent of the renormalization scale p. Since we know the renormal- 
ization group equation for q, q and the soft factor, we can easily derive the renormalization 
group equation for the hard part, 



dH(Q 2 /p 2 ,p) \\TT/rfl/ 2 \ 

^ dp = ~ ^ lF ~ 1S ^ P '' ^ IV ,P) ■ 



(76) 



The solution is 



H(Q 2 /p 2 ,p) = exp|-^^[4 7F - 75 (p)]}iJ(g 2 //i' 2 ) p) • (77) 

To factor out the large renormalization logarithms, one can choose p to be at low scale such 
as pl, and p! at high scale such as pn- Therefore, we write, 

H{Q 2 /plp) = exp (- ^[4 7F - 75 ] j H(Q 2 /C 2 x(,p) , (78) 

where H(Q 2 /C 2 x(,p) contains no large logarithms. 

Collecting the above results, one has both the renormalization and soft-collinear loga- 
rithms summed in the following expression: 

F(x B ,z h ,b,Q 2 ) = q (x B ,z h b,p 2 L ,p L /C 2 ,p) q (z h ,b, p 2 L , p L /C 2 , p) S(b,p 2 L ,p)H (l/C 2 p,p) 



x exp < -2 



c ^d_p 



P 



1k{ol(p)) -K(b,p L ) 



-G(p/C 2 ,p)-2 lF + - ls (p) 



(79) 
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where all large logarithms have been factorized in the exponential factor. For the physics 
discussed in this paper, we do not want significant large logarithms, because otherwise the 
transverse momentum of the hadron yield is generated mostly by soft-gluon radiations. To 
avoid them, Q 2 can only be moderately large compared to Ph±- On the other hand, in 
this kinematic regime, the contributions from power-suppressed terms might not be entirely 
negligible. 

Finally, the choice of p. According to its definition, we must have p > 1 although the 
physics is independent of p. However, if p is too large, one has large logarithms in the hard 
part and the convergence of the perturbation series might be spoiled. Therefore in practice 
one might choose a p, for example, somewhere in between 3 and 10. 



VI. CONCLUSION 

In this paper, we showed that a factorization theorem exists for semi-inclusive deep- 
inelastic scattering with detected hadron momentum P±h «C Q. Ph± can either be soft, i.e., 
on the order of Aqcd, or in the perturbative domain Aqcd- We have mainly focused 
on the former case although the result is valid also for the latter. For P±h ^ Aqcd, the 
theorem can be simplified by an additional factorization of the TMD parton distributions 
and fragmentation functions 0. 

We established the theorem by first considering the example at the one-loop level. In this 
case, the calculation of the parton distribution, fragmentation function, soft factor, and the 
SIDIS cross section was straightforward. The example demonstrated that the factorization 
indeed works. 

At higher-order in perturbation theory, one can use the formalism developed by Collins, 
Sterman, and Soper and others. Starting from the most general reduced diagrams, we 
showed the factorization of collinear gluons from the hard part and the soft gluons from the 
collinear part, matching the jet factors with the distribution and fragmentation functions. 
The factorization scale p and the p independence of the physical cross section allows one to 
sum over large logarithms involving scales Q and Ph±- 

The present result can be easily extended to the situations where the target is polarized 
or the polarization of the final state hadron is measured. It can also be extended to the 
case where the transverse momentum of the hadrons are integrated with a weighting factor. 
Finally, all results here can be obtained also in the framework of the soft-collinear effective 
theory. 
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FG02-93ER-40762. J.P.M. was supported by National Natural Science Foundation of P.R. 
China (NSFC). X. J. is also supported by an Overseas Outstanding Young Chinese Scientist 
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